Electric circuit networks equivalent to chaotic quantum billiards 
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We formulate two types of electric RLC resonance network equivalent to quantum billiards. In 
the network of inductors grounded by capacitors squared resonant frequencies are eigenvalues of the 
quantum billiard. In the network of capacitors grounded by inductors squared resonant frequencies 
are given by inverse eigen values of the billiard. In both cases local voltages play role of the wave 
function of the quantum billiard. However as different from quantum billiards there is a heat power 
because of resistance of the inductors. In the equivalent chaotic billiards we derive the distribution 
of the heat power which well describes numerical statistics. 



PACS numbers: 03.65.Ge,03.65.Yz 



I. INTRODUCTION 



Firstly electric circuit models representing a quantum particle in the one-dimensional potential 

h 2 d 2 Tp{x) 



2m dx 2 



■V(x)ip(x) =Eip(x) (1) 



were considered by Kron in 1945 Three types of equivalent circuits were established. The first one contains 
positive and negative resistors and in each state the currents and voltages are constant in time. The second and third 
models are similar and consist of inductors and capacitors and the currents and voltages are sinusoidal in time. Here 
we consider the stationary Schrodinger equation in two-dimensional billiards in hard wall approximation 

- V 2 ^(x,y) =a/ ) (x,y), (2) 

where the Dirchlet boundary condition is implied at the boundary C of the billiard: 

Mo = 0- (3) 

Here we use Cartesian coordinates x, y which arc dimcnsionlcss via a characteristic size of the billiard L, and corre- 

spondingly 6=^,^0 = a£*. 

There is a complete equivalence of the two-dimensional Schrodinger equation for a particle in a hard wall box to 
microwave billiards 0] . A wave function is exactly corresponds to the the electric field component of the TM mode of 
electromagnetic field: if)(x, y) <-> E z (x, y) with the same Dirichlet boundary conditions. This equivalence is turned out 
very fruitful and allowed to test a mass of predictions found in the quantum mechanics of billiards Q • On the other 
hand, models for the equivalent RLC circuit of a resonant microwave cavity exist which establish the analogy near an 
eigenfrequency Q • Manolache and Sandu Q proposed a model of resonant cavity associated to an equivalent circuit 
consisting of an infinite set of coupled RLC oscillators. Therefore, there to be a bridge between quantum billiards 
and the set of coupled RLC oscillators 5]. In fact, we show here that at least, two models of electric resonance 
circuits (ERC) can be proposed. In the first model shown in Fig. Q]the eigen wave functions correspond correspond 
to voltages and eigen energies do to squared eigen frequencies of ERC. In the second model shown in Fig. [3 the eigen 
energies of quantum billiard correspond to the inverse squared eigen frequencies of the electric network. The electric 
network analogue systems allow to measure not only typically quantum variables such as probability and probability 
current distributions but also a distribution of heat power in chaotic billiards. Moreover intrinsic resistances of the 
RLC circuit allow to model the processes of decoherence. 



II. ELECTRIC RESONANCE CIRCUITS EQUIVALENT TO QUANTUM BILLIARDS 



If to map the two-dimensional Schrodinger equation onto numerical grid (x,y) — ao(i,j),i — 1,2, ...N x ,j = 
1,2,... N y one can easily obtain equation in finite clement approximation 

+ i>i,j-i + ipi+ij + ipi-i,j + (<4 E - 4)V^.j = 0. (4) 
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The equivalent Hamiltonian is the tight-binding one 



# = -EEiy'Xitf +b i 



(5) 



where vector b, |b| = 1 runs over the nearest neighbors. 

Let us consider the electric resonance circuit shown in Fig. ^ Each link of the two-dimensional network is given 
by the inductor L with the impedance 



Zl = i^L + R 



(6) 



where R is the resistance of the inductor and u is the frequency. Each site of the network is grounded via the capacitor 
C with the impedance 



1 



zc 



(7) 




FIG. 1: The first model of resonance RLC circuits. 
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FIG. 2: The second model of resonance RLC circuits. 



The Kirchoff 's current law at each site of the network gives 

— [Vi. ]+ i - Vi tj + Vij-! - Vij + Vi+ij - Vij + Vi-ij - Vij] - — Vij = 0, 
zl z c 



(8) 



where Vij are values of voltage at the site One can see that this equation coincides with the discretized version 

of the Schrodinger equation with the eigenenergies as 



<yJ' 2 1 

zc 



aik 2 = -— = LClo 2 - iRCuj = — i 



— % 2", 

^0 



(9) 
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where ujq = l/VLC and 7 = R/L are the eigen frequency and the linewidth of each resonance circuit. 
For the second network of electric resonance circuits shown in Fig. [5] we obtain 

— [Vij+i - Vi,j + Vi,j-i - Vi. 3 + V i+ t,j - Vi,j + Vf-1,3 - V itj ] - —Vi,j = 0. (10) 

ZC Zl 

Therefore, comparing with we have 

^ = -? = W + £ = ^+^- (ID 

where 7 = RC . This network is interesting in that its eigen frequencies are inverse to the eigenenergies of the quantum 
billiard. 

There are many ways to define the boundary conditions (BC). Let it be (iB,js) are sites which belong to the 
boundary of the network. If these sites are grounded, we obtain obviously the Dirichlet BC © V\ B = 0. It they are 
shunted through capacitors we obtain the free BC (the von Neumann BC). At last, if the boundary sites are shunted 
through resistive inductors, the BC correspond to mixed BC. 



III. ANALOG OF THE CHAOTIC BUNIMOVICH BILLIARD 



A real electric circuit network has three features which can make a difference if to compare to the quantum billiards. 
These are 1) a discreteness of resonance circuits, 2) tolerance of electric elements, and 3) resistance of inductors. In 
practice the discreteness has no effect for A > 10ao where A is a characteristic wavelength of wave function, and ciq 
is the elementary unit of the network. Numerically we consider the electric network with shape as a quarter of the 
Bunimovich billiard. The distribution of real part of the wave function of the billiard mapped on the electric circuit 
network with ciq = 1/100 is shown In Fig. |3| (a). The wavelength A = 27raoujo/u} = 0.115 for parameters given in 
caption of Fig. |3| We take the width of the billiard as unit. One can see distinctive deviation from the Gaussian 
distribution which is result of multiple interference on discrete elements of the network. 

It is known that a noise, for example, temperature, smoothes the conduction fluctuations for transmission through 
quantum billiards 0, 0|- In present case the tolerance of circuit elements, capacity and inductance, plays role of 
the noise. Therefore we can expect that increasing of the tolerance can suppress fluctuations of the distribution of 
the wave function of the discrete electric circuit network. In fact, even the 1% tolerance substantially smoothes the 
distribution of the wave function as shown in Fig. [3] (b)-(d). We consider that the fluctuations of capacitors and 
inductors are not correlated at different sites. 

Finally we consider as a damping caused by resistance of the inductors effects the distribution of the wave function 
in the electric RLC resonance circuit network. In order to excite the network we apply external ac current at single 
site of the network. Fig. 0] shows the probability density in the quarter of the Bunimovich billiard for two values of 
the resistance R. One can see from Fig. 0J (right) a localization effect because of a damping of the probability density 
flowing from ac source (see also Fig. 0). The characteristic length of space damping can be easily estimated from Eq. 
© which gives us 

The distributions of the probability density p — \V\ 2 for open quantum chaotic billiards were considered in many 
articles H, 0, 0, 0, 0] for the case of zero damping. Here we follow 0, Q] and perform the phase transformation 
V — > V exp(i#) = p + iq which makes the real and imaginary parts of the wave function V independent. Introducing 
a parameter of openness of the billiard |14j 

e 2 = f 2 . (13) 

where a 2 = (p 2 ), o 2 q = (q 2 ) we can write the distribution of probability density as ^3 

f(p) = Mexpf-MVKoO^/ 5 )) (I 4 ) 

with the following notations 

Hfi +e V=UH' (i5) 
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u/(u 2 > 1/2 u/(u 2 ) 1/2 

FIG. 3: The distribution of real part of the wave function of the quarter Bunimovich billiard mapped on resonance RLC circuit 
with elementary unit ao — 0.01, uj = 1.722MHz, L — 0.1 mH, C = 1 nF, R — 0. (a) There is no tolerance of the electric circuit 
elements, (b) The tolerance equals to 1%. (c) The tolerance equals to 3%. (d) The tolerance equals to 5%. Each distribution 
in (b) - (d) is averaged over 100 realizations of the electric network. 

and Iq(x) is the modified Bessel function of zeroth order, This distribution is shown in Fig. EJby solid lines while the 
Rayleigh distribution f(p) — exp(— p) is shown by dashed lines. The Rayleigh distribution specifies the distribution 
of completely open system. One can see from Fig. (a, b) that the statistics of the probability density follows the 
distribution l|14|) irrespective to resistance R. However with growth of the resistance the distribution i|14[) tends to 
the Rayleigh distribution (Fig. [S](c, d)). Since the larger resistance the more quantum system is open, this tendency 
of statistics of the probability density is clear. 




FIG. 4: Views of probability density of the quarter Bunimovich billiard mapped on resonance RLC circuit with elementary 
unit a = 0.005, w = 0.8611 MHz, L — 0.1 mH, C = 1 nF. Left R = 0.5 fi, right R — The point of connection of external 
ac current is at maximum of the probability density. 
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FIG. 5: (Color online) Distribution of probability density of the quarter Bunimovich billiard mapped on resonance RLC circuit 
with the same parameters as given in Fig. |U (a) R = 0.1 f2, Q = 3162, e = 0.2488, (b) R = 0.3 fi, Q = 1054, e = 0.5308, (c) 
R = 0.5 Q, Q = 632, e = 0.6996 and (d)iJ=lfi,Q = 316, e = 0.9164. The distribution (TQ is shown by solid red line, the 
Rayleigh distribution f(p) = exp(— p) is shown by dashed green line. 
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(a) Re(l x )/(Re(l x ) 2 > 1/2 (b) P / (P > 

FIG. 6: (Color online) (a) Statistics of the real part of the x-component of electric current I x compared to the Gaussian 
distribution shown by solid red line, (b) Statistics of the heat power compared to the distribution 1281 shown by solid red line. 
Here the quarter Bunimovich billiard is taken with u> = 1.163 MHz, R = 0.1 £1, L = 0.1 mH, G = 1 nF. 



IV. THE HEAT POWER 



In open systems the probability current density corresponds to the Poynting vector. The last equivalence allowed 
to test in particular universal current statistics in chaotic billiards [T^. H3 | . However in the electric resonance circuit 
there is the heat losses because of the resistance. A local power of the heat losses is denned by formula 0] 

P = §[Re(Ix) 2 + IHIx) 2 + Re(Iy) 2 + Im(Iyf] = f [%? + \Iy\% (16) 
where I x , I y are local components of the electric power flowing between sites of the electric network: 

RI x (i,j) = V i+ld -Vi,j, RI y (i,j) = Vi, j+l -Vi,j. (17) 
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If to approximate the true state with the Berry conjecture 

v(x, v) = ^2 a j ex p[* ( k j r 



(18) 



where aj and <pj are independent random real variables and kj are randomly oriented wave vectors of equal length, 
then V is the complex random Gaussian field (RGF) in the chaotic Bunimovich billiard. The derivatives of V are also 
independent complex RGFs. The components I x , I y form two complex RGFs with the probability density of these 
fields 



111 _|_ 77/2 



I I 1 Mi + 1 4 

f{l' x j' y X,Iy) = ^2^2 eX P j~2 y X ~2 " + - -2 " J j 

where £ = fle(7 x ), = i?e(U = M{»), I" = Re(I v ), a? = (I' X 2 ),(I?), ^ = (If 
computations we use that average over the billiard area 

is equivalent to average over three complex GRFs 

d 2 Vd 2 I x d 2 I y f(Re(V),Im(V))f(I' x ,I'I x \ /").... 



(19) 

In numerical 
(20) 

(21) 

An example of the distribution of the real part of I x is presented in Fig. El (a) which shows that numerically this 
value is, in fact, the RGF. A definition of the probability distribution (|19|) is relied on that the Berry function 
(|18fl is isotropic in space : (I' 2 ) = (I y 2 ), (I x 2 ) = (Iy 2 )- I n fact, an anisotropy of the shape of billiard effects an 
anisotropy. However this effect is the boundary condition's one which has of order LpX/A ~ A, where Lp is a length 
of the billiard perimeter, and A is a characteristic wave length of wave function in terms of the width of the billiard. 
Therefore, for the excitation of the eigenfunction with sufficiently high frequency we can use the distribution function 
(|19H . Table 1 of numerically computed mean values confirms this conclusion. 

Table 1. Numerically computed mean values. 



u, MHz 


the wavelength A in terms of the billiard's width 






e 






0.8611 


0.1154 


0.095 


-0.128 


0.2488 


1.1623 


0.0854 


0.056 


0.050 


0.6103 



To find the distribution of the heat power 116JI it is convenient to begin with a characteristic function 
9(a) = (exp(* a P)) = J d 2 I x d 2 I y f(r x ,r y ,r x \I^)eM™R{\I x \ 2 + \I y \ 2 }/2). 
Substituting (|T^|l we obtain 

ew = _(fi±f2! ! 

" -W (. + ,^t)(. + ,^a)- 

A knowledge of the characteristic function allows to find the heat power distribution function 

1 f°° 2/i 

— / dae(a)e X p(-iaP) = ^-exp(-vP/(P))smh(vP/(P)), 
2tt J.oo u{P) 



(22) 



(23) 



f(P) 



where formulas l|15l take the following form 



2c?af ' 



2o*o? 



(24) 



(25) 
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For a 2 w of the distribution takes the very simple form 

/( P ) = tJ^ cx p(- 2F /< F ))' ( 26 ) 

Even for this case the distribution of heat power differs from the distribution of the probability current 14]. The 
parameter of openness of the billiard can be approximated as 

6 2 = 4- (27) 
of 

It is easy to obtain from the Schrodinger equation that 2 of = Ea 2 , 2 of = Ea 2 from which the last equality in 127|) 
follows. 

Then the heat power distribution function (|24|l can be written as follows 

, s 1 + e 2 f / (l + e 2 )P\ / (l + e 2 )P\) , s 

This distribution is shown in Fig. (b) which as one can see nicely describes numerically computed statistics of the 
heat power. If to introduce a value 

((P-(P)) 2 ) 

P = Jpy ' (29) 

then one can derive the relation between this parameter and the parameter of openness i|27|) 

(30> 

If the quantum system is fully opened, e = 1, and we have from l|3U|) that a 2 p = 1/2. For the limit of closed quatum 
system we obtain correspondingly that <7p = 1 . 



V. SUMMARY AND CONCLUSIONS 



We established two types of the electric circuit networks of the RLC resonant oscillators in which voltages play role 
of quantum wave function. Specifically we considered the electric networks with the Dirichlet boundary conditions 
which are equivalent to the quarter of the Bunimovich quantum billiard. In fact, the electric circuit network has three 
features which can make a difference if to compare with the quantum billiards. These are a discreteness of resonance 
circuits, tolerance of electric elements, and resistance of inductors. We showed numerically that the first two features 
are conceal each other. The resistance of the electric network gives rise to a heat which can be described locally by 
the heat currents. Assuming that the wave function in the billiard can be given as the complex random Gaussian 
field we derived the distribution of the heat power which well describe numerical statistics. 

The third feature of the electric network, resistance has principal importance. The resistance of the electric network 
is originated from inelastic interactions of electrons with phonons and other electrons which give rise to irreversible 
processes of decoherence. With growth of the resistance the wave function is becoming localized. We studied as the 
probability density and the probability currents evolves with increasing of the resistance. Therefore we can conclude 
that the resistance violates the equation Vj = 0. In fact Fig. [7]demonstrates unusual behavior of quantum streamlines 
[rH with growth of the resistance. One can see that the quantum streamline are terminating at the vortex cores. 
The vortices serve as sinks for the probability density shown in Fig. [7| (top) as spots. Therefore the resistance of the 
inductors in the equivalent electric networks is simple mechanism of a deterioration of ballistic transport similar to 
the Biittiker mechanism [Isj . 
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FIG. 7: Top: quantum streamlines in the quarter of the Bunimovich billiard flowing from point shown by star at which the 
external ac current is applied. (Color online) Bottom: zoomed part of top figure. Blue lines show the streamline, red and green 
lines are the nodal lines of the real and imaginary parts of wave function correspondingly. The points at which the nodal lines 
intersecting are centers of the vortices |l7|. The wave function corresponds to Fig. El (right) with the same parameters. 
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